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Optimal quadrature formulas of closed type in the space L2 (0, 1) 
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Abstract 

It is discussed the problem on construction of optimal quadrature formulas in the sense of Sard in the 
space Lj (0, 1), when the nodes of quadrature formulas are equally spaced. Here the representations of 
optimal coefficients for any natural numbers m and A^ are found. 
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< 

c^ ■ Key words and phrases: Sobolev space, optimal quadrature formula, positive coefficients, 

error functional. 

T-H ■ 
>, 

f^ . Consider quadrature formulas of the form 

1 N 



^{x)dx = Y,C[P]^[P] 



O : ;i^ /3=o 

o 

with the error functional 

N 

>< : (i{x) = qo,i](x) -J^Cim^ - h^), (2) 

C^ '_ 13=0 

where C[/3] are the coefficients of formula (1), [/3] = /i/3, h = 1/N, N = 1,2, 3, ..., the function 
V9(x) belongs to the space L2 (0, 1). Norm of functions in the space L2 (0, 1) is determined 
by formula 




||v.(x)|4") (0,1)11= y {y^^'-\x)) dx 
Below [13] = h(5. 

^This article is the fifth section of the author's Candidate dissertation entitled "Optimal formulas of approximate integration 
for differentiable functions" - Novosibirsk, 1983. 140 p. 



Quadrature formulas of the form (1) is closed type, since the points a; = and x = 1 are the 

nodes of the formula. 

The difference 

} 

{i{x),v{x)) = / ^{x)dx-j2cmvm 

/5=o 

is called the error of formula (1). 

The problem of construction of optimal quadrature formulas (1) in the sense of Sard in the 

space L2 (0, 1) consists of computation the quantity 

||IW|4""-(0,l)ll = inf sup l'^'"'f'"''l . (3) 

where L^ (0, 1) is the conjugate space to the space L2 (0, 1). 

The coefficients C[(3] satisfying equality (3) (if there exist) are called optimal and are denoted 

o o 

as C [/?]• Formulas of the form (1) with coefficients C [/?] are called optimal quadrature formulas 
in the sense of Sard. 

The problem of construction of optimal quadrature formulas in the sense of Sard in different 
spaces are investigated by many mathematicians. In the space Lg this problem was inves- 

o 

tigated by S.L.Sobolev in [1] and for optimal coefficients C [/?] the system of linear equations 



was obtained, which for X[s = [/3] = h/3 {/3 = 0, N) has the form 

l^^^^U.{2m-l)\ +Z^^"[^] = 2i2m)\ ' [^] ^ f^' ^1' ^4) 

7=0 ^ ^ 0=0 ^ ^ 

J2cb]br = ^^ a = 0,m-l. (5) 

o 

In this system C [/?] and Aq,, a = 0, 1, ..., m — 1 are unknowns. 

Note that the existence and uniqueness of the solution of the system (4)-(5) was proved by 
S.L.Sobolev in [1]. 

o 

The aim of this work is finding the explicit forms of the optimal coefficients C [/?] • 

o 

Later for convenience optimal coefficients C [P] we remain as C[/3]. 



Here mainly is used the concept of discrete argument functions and operations on them. For 
completeness we give some definitions about functions of discrete argument. 

Assume that v^(x) and iplx) are real- valued functions of real variable and are defined in real 
line R. 

Definition 1. Function ^{hp) is called the function of discrete argument^ if it is given on 
some set of integer values of /3. 

Definition 2. By inner product of two discrete functions ip{hP) and ip{hp) is called the 
number 

oo 
/3=-oo 

if the series on right hand side of the last equality converges absolutely. 

Definition 3. By convolution of two functions ^{h(5) and ip{h(3) is called the inner product 

CO 

7=— oo 

Suppose C[/3] = for /3 < and (3 > N . Then, using definition 3, system (4)- (5) can be 
rewritten in the following form 

Gn^AP] * C[(3] + P„_i[/3] = /„[/?] when [f3] e [0, 1], (6) 

C[/3] = when [/3]^ [0,1], (7) 

^ 1 

VC[/3] [/?]" = ^, a = 0,l,2,...,m-l, (8) 

where Pm-i[/3] is polynomial of degree m — 1 with respect to [/3] and 

_ [/3]^-^sign[/3] 
^-■^f^J- 2(2m-l)! ' ^^^ 

^'"^^J 2(2m)! ■ ^^' 

Consider the following problem. 

Problem A. Find the discrete function C[f3] and unknown polynomial Pm-i[P]- 

Denote 

v[/3] = G^M*C[/3] (11) 



and 



u[l5]=v[P]+P^^r[P]. (12) 

Then we have to express C[/3] with the help of «[/?]. For this we must construct discrete operator 
Dm[/3]j which satisfies the equation 



hD^[l3]*G^M=m. 



(13) 



where (5[/3] equals when /3 7^ and equals 1 when /3 = 0, i.e. 5[[5] is the discrete delta-function, 
Cm.ii/S] is defined by formula (9). 

In connection with this in the work [2] the discrete analogue Dm[(S\ of the differential operator 
(P™- /dx^^ was constructed. 

In [2] the following was proved. 

Theorem 1. The discrete analogue of the differential operator d'^™' / dx^"^ have the form 



Dr. 



(2m- 1)! 



E 



;i - Qkf 



Qk 



1 QkE- 



2m--l\(ik, 



m— 1 



l + Z 



^-qk 



v2m+l 



^ E2m~l{<ik) 
,2m-l , Y^ U Ik 



for 
for 



>2, 



(14) 



-^ + E 



i2m+l 



fc=l 



qtE; 



for /3 = 0, 



2m-l\qk 



where E2m--i{x) is the Euler-Frobenius polynomial of degree 2m— 1, qk are roots of the polynomial 
E2m~2{x), \qk\ < 1; h is small parameter. 

Theorem 2. For the operator -Dm[/3] and monomials [/3]'' = {hPY the following are true 



/3 



when < k < 2m — 1, 

(2m)! when k = 2m, 



(15) 






/i^"(4m)!E2m 
(2m)! 



when 2m + 1 < k < Am — 1, 
when k = 4m. 



(16) 



Taking into account (13) and theorems 1, 2 for optimal coefficients we have 



C[/3] = hDm[/3]*u[/3]. 



(17) 



Thus, if we find the function u[l3], then optimal coefficients will be found from (17). 

In order to calculate convolution (17) we need representation of the function u[P] in all 
integer values of (3. From (6) we have u[f3] = /m[/3] when [/3] G [0, 1]. Now we must find 
representation of ■u[/3] for /3 < and P > N. 

Since C[/3] = when [(3] ^ [0, 1], then 

C[f3] = hD^[f3]*u[f3] = 0, [/3]^[0,1]. 

We calculate the convolution v[l3] = Gm,i[l3] * C[l3] when [/3] ^ [0, 1]. 
Suppose /3 < 0, then taking into account (4), (5) we have 

2^(j + l)!(2m-l-j)! 2^j\{2m-l-j)\2-. ^^^^^^ ' ^ > 

Denote 

2 ^ (j + l)!(2m - 1 - j)! 2 ^^ j!(2m - 1 - j)! ^ 

Then from (18) for v[f3] when /3 < we obtain 

t'[/3] = -g2m-i[/3]-i?™-i[/3]. (19) 

Similarly for the case P > N we have 

v[P]=Q2m^im + Rm-im- (20) 

Denoting 

i?;;-i[/3] = ^r„-i[/3]-^m-i[/3], (21) 

R+_,m = Pm-im+Rm-im, m 

and taking into account (19), (20), (12) we get following problem. 
Problem B. Find solution of the equation 

hD,^[/3]*u[/3] = 0, [/3]^[0,1] (23) 



having the form 

' -Q2„.-i[/3]+i?-_i[/3], /3<0, 

«[/?]= \ /™[/3], 0</3<iV-l, 

g2™-i[/3]+i?+_i[/3], /3>iV, 
v 

where R^^_i[(3] and i?^_]^[/3] unknown polynomials of degree m — 1. 

If we find i?;;_i[/3] and i?+_i[/3], then from (21), (22) we get 

„ _, 1 
2 
1 
2 



(24) 



m-l 



i?™-i[/3] 






Unknowns R 



m—l [ 



and -R^_]^[/3] can be found from (23), using the function Dm[f3]. Then we 



obtain exphcit form of ■u[/3] and from (17) will be found the optimal coefficients C[(3]. Thus, 
the problem B and respectively the problem A will be solved. 

But here we will not find i?~_j^[/3] and R'^_i[(3], instead of them, using -Dm[/3] and taking 
into account (17), we will find expressions for the optimal coefficients C[/3]. 

The main result of the present work is the following. 

Theorem 3. Optimal coefficients of quadrature formulas of the form (1) on the space 
L2 (R) have the form 



m—l 



2 /L^ K i_q^ 



C[(5] = h < 



k=l 
m—l 



i+UdkU + it' 

k=l ^ 



for I3 = 0,N. 

for 13 
for/3<0,/3>N, 



l,2,...,iV-l, 



where 



m—l 



E*E 



Qk 



.ly+lqN+^ 



(qk - iy+' 



A^O^' 



B 



i+i 



J + 1- 



J 



1,2,..., m-l. 



k=l i=l 

Qk are roots of the Euler-Frobenius polynomial E2m-2{x) of degree 2m — 2. 
In the proof of theorem 3 the following preliminary results are used. 
The following formula is true [3] 

n—l -, k / \ i „ k 

q 



7=0 



q^[l 

^ 1=0 ^ 



q 



A^O'^' - 



q" 



q 



7=0 ^ 



q 



q 



AS''! 



(25) 



where A*7'^ is finite difference of order i from 7", A'O'^ = A*7"|^=o- 

And also we use the following well-known formula (see, for example, [4]) 

/3-1 fc+l , , 



It is known in [5], that the Euler-Frobenius polynomials Ek{x) have the form 

xEkix) = (1 - xf-^^D'-^, (27) 

1 — xr 



where 

D = x^, D' = x^D'-\ 

ax ax 

In [5] it was shown that all roots g- of the Euler-Frobenius polynomials Ei^{x) are real, negative 
and distinct: 

gf ) < gf < ... < gf < 0. (28) 

Furthermore for the roots (28) the following is true: 



1. 



1j ■ ?fc+l-j 
k 



If we denote Ek{x) = "^ as x*, then coefficients as of the Euler-Frobenius polynomial Ek{x) 

s=0 

are expressed by formula 

This formula was obtained by Euler. 

From definition of Ek{x) we get following statements. 

Lemma 1. For the polynomials Ek{x) the following recurrence relation holds 

Ek{x) = {kx + l)Ek-i{x) + x(l - x)El_^{x), (29) 

where Eo{x) = 1, k = 1,2, .... 

Lemma 2. The polynomial Ek{x) satisfies the identity 

Et,{x)=x''EJ^ (30) 



or otherwise ai = %_,,, s = 0,1,2, ...,k. 
Lemma 3. Polynomials 







Pk{x) 


= 


{x~ 


fe+i 

i=l 


(x - 1)^ 


> 1 


{x 


/ \ 


'1 

X 


-1 


k+1 fc+1 . 


X Y 


k 1 


1 — X J 



(31) 

^ IX- - Lr 

and 

/l\ /I \ fe+i fc+i / ^ \i 

A*0^'+^ (32) 

are i(/ie Euler-Frobenius polynomials Ek{x) and E^ (-) respectively. 
Lemma 4. T/ie operator D.fn[[i] satisfies the relation 

oo 

[D^m, W^] = E ^™[/5][/5]'" = (2^)!- (33) 

/3=-oo 

Proofs of Lemmas. 

Proof of lemma 1. From (27) obviously that 

Eu.,{x) = x-\l - xf+'D^-^—^—. (34) 

(1 — xY 

Differentiating by x the polynomial Ek_i{x), we obtain 

k^~2n„^ I 1 \ nfe-i ^ I ^k[X) 



E',_,{x) = -{l-x)'x~\kx + l)D 

(1 — X)'' X[l — X) 

Hence and from (34) implies that 

{kx + l)Ek-i{x) + x{l - x)E'k_-^{x) = {kx + l)a;-^(l - x)''+^D''-^ 

-(1 - x)''+'x-\kx + l)D'-'—^^ + Ek{x) = Ek{x). 

1 — xr 



X 



{i-xy 



Thus lemma 1 is proved. 

Proof of lemma 2. Lemma 2 we will proof by induction. For k = 1 from (27) we find 

Ei{x) =x + l. 

Assume that for /c > 1 the equality an~ = a^Siln^ n = 0,1, ...,k — 1 is fulfilled. Suppose that 
On = for n < and n > k — 1. 



From (29) we get 



af = {s + l)af-^) + (A; - s + l)af_"/^ 



then, using assumptions of the induction, we obtain 

4'i = {k-s + l)at/) + (. + l)at^'l, = {k-s + l)af_-/) + (. + l)af-^) = af , 

and Lemma 2 is proved. 

Proof of lemma 3. Consider relations (29) 

Ek{x) = {kx + l)Ek^i{x) + x{l - a;)_E^_;^(a;), 

Eo(x) = 1. 

If Pk{x) also satisfies this relations then lemma 3 will be proved. From (31) evidently that 
Po{x) = 1. We denote 

Vk{x) = {kx + l)Pk^i{x) + x{l - x)Pl^^^{x). 

Using the equalities 



i=l 



PUix) = k{x - If-' V j^^ -{x- If V -^^^, 
after some simplifications we have 

Vfc(x) = {kx + l)Pk-i{x) + x{l - x)Pl_^{x) = 

= Pu-i{x) + x{x - If J2 -r-TTY ^ ^ ^'°'^^ + "^'^^'^ ~ ^^'"- 
Doing change of variables a; — 1 = y, we obtain 

k 

Vk{y + 1) = Y,y''\^'^' + A*-''0^)(^ + 1), 

i=0 

k 



Pfe(l/ + l) = 5^/-W+i0'=+^ 



1=0 



Consider coefficients of the polynomial Vk{y + 1) 

i 



a=0 



E(-i)- :h^+E(-i 



j+i 



\i+l— a 



a=0 



i + 1 



a 



a 



(^ + 1) 



a=0 



a=0 ^ 



(i+1) 



5:(-i: 



a=0 



a 



a=0 



+ Z(-i)'*'"°l I " l + l " . I I "' ■ + (' + 1) 



a/ \a — \ 



t^ + ll 



a=0 



$:(-ir-"( ^Ja^ + (. + l 



(^ + 1) 



i+1 



E(-i)^^''" 



a=0 



^ + 1 ^ fc+1 _ ^i+lQfc+1 



a 



a 



Hence we get Vk{y + 1) = Pfc(t/ + 1), i.e. 



Pk(x) = EJx) 



Lemma 3 is proved. 

Proof of lemma 4. Using (14) we have 



oo m— 1 



(oo m— i /-. \2m+l 



+2[1] 



2'"') =2-(2m-l)! (^ 



QkE2m-l{Qk) 



Hence, by virtue of well-known formula ^ q"''-)^ = j^- X]j=o ( T^ ) ^^O'^) we obtain 

7=0 '^ V 'J/ 



''m—l 



2m+l 



According to lemma 3 



1 ^"^ / 
1 - at ^ I 1 



feJ -^ ^ ^fc ~o ^ ^ ^'^ 



Qk 



^^QZm ^ ^ 



'm— 1 



[DM, Wn = 2(2m - 1)! ( 5^ '^ : ^'^'!""'! ■ .. ^''Tr^^^^^'l...,. + 1 



^ qkE2m-iiqk) (1 - qk){Mqk - 1) 



2m 



10 



Now, using lemma 2, finally we get 

[D„M. \B?"\ = 2(2,„ - 1)! fg ''-""'""i ■ 71 "i^^^'f ,„„. + 1 

Vfc^ ^kE2m-i[qk) (1 - gfejgfc (1/gfc - 1)^™ 

= 2-{2m-l)\-m = (2m)!. 
Lemma 4 is proved. 

Proof of theorem 3. 



Suppose (3 = 1,N — 1. Then from (17), using (24) and definition of convolution of discrete 
functions, we have 

oo ^ —1 

C[P] = hD^[P]*u[P] = h Y, D^[P-^]uYi] = hl Yl Dm[/3-l]{-Q2m-lb]+Rn.-lb]) + 

7=— oo ^ 7=— oo 

N oo 

+ J2DrrAP-l]fmb]+ Y. Dm^ - l] {Q2m-lb] + K^Al]) 
7=0 7=Ar+l 

(oo oo 

Y Dm[P - l]frnb] + ^^-t^^ + 7] (-^2™-! [-7] + i?™-l[-7] " /m[-7]) + 
7=— 00 7=1 

00 \ 

+ 5^ D™[iv + 7 - /3](g2™-i[iv + 7] + i?™-i[iv + 7] - ^[iv + 7]) I 
7=1 ^ 

(00 
D^[I5] * Um + J]l^™[/3 + 7](-Q2„^-i[-7] + Kn^A-l] - U-l]) + 
7=1 

00 \ 

+ 5^Z}™[iV + 7-/3](Q2™-i[iV + 7] + i?+-i[iV + 7]-/™[iV + 7])). 
7=1 ^ 



7 
Hence, using (14), (15) and (33), we get 






K=l 7=1 

m—1 



/o -|\i /I N2m+1 °" \ 

+ E^r' ,2: ,r^\,, E^n^^»--i[^+7] + ^;t.-i[iv+7]-/4iv+7]) . 

^^ Ai qkJ^2m-l[Qk) ~[ J 

(35) 



11 



We denote 



oo 



^Z' (36) 

7=1 



where k = l,m — 1. Then when /3 = 1, A^ — 1 for the optimal coefficients 



m— 1 



C[P]=h\l + Y, (duqi + p.gf -'^) j . (37) 

Now from (8) when a = 0, 1 keeping in mind (37) for C[0] and C[N] we obtain the expressions 



C\N] = h(- + Y^ (MLZJ^]^ + Kq^-^^ - qk){dk - Pk) \ \ 



( 1^2 //' (38) 

[qk-iy J) 

(39) 



where dk and pk are determined from (36). 



Now in order to prove theorem 3 sufficiently to show the equalities dk = Pk, k = l,m — 1. 



Now we will prove, that dk = Pk, k = l,m — 1. 

Consider the convolution in equality (6) and we rewrite it in the form 



9[^]=G^A^]*C[^] = J2C[l] 



where 



2(2m-l)! 

7=0 ^ ^ 



First we consider Ti. For Ti using (34) we have 



,lS^^^(-S(«--"i)^^^fe^ 



7=1 \ k=l 

^f°](2;;7^i)! + (2;;7^i)!(,E^ + 



7=0 

12 



m—1 



/3-1 



13-1 



A;=l V 7=0 7=0 / ^ 

Using by formulas (25), (26) the expression for Ti we reduce to the form 



Ti = C[0] 



(2m -1)! (2m -1)! 



^ {2m-l)\B2„,-j ^ 
^ J!(2m-j)! ^ 



m— 1 



+ E ^^^^ 



k=l 



Qk 



- 1 ^ fOA. - D* Qk-l ^ {Qk - lY 



i=0 



+Pkqk 



2m-l 



1-?^ ?? Vi 



gfe 



gfc 



A*o 



in2m- 



gfc 

/3 2m-l 



i=0 



+ 






Qk 



Qk 



A'f3 



i o2m~l 



(42) 



i=0 \ -"-/ -">- j=o 

Keeping in mind that q^ is the root of the Euler-Frobenius polynomial E2m-2{x) of degree 
2m — 2, using lemma 3 and equality (38), after simplifications, from (42) we get 



Ti 



h 



•2m 



(2m -1)! 



/3 



2m 



2m 



+ 



2m- 



2m-2 



m—1 2m— 1 , I / 1 \j N+i 

S^ yz,IIL- l)'.lD2m-j ai ,\^\^ -dkQk + l-i j Qk Pk ^i p2m-l 



(2m-l)!i?2m-, ,, 

^j^ j!(2m-j; 
Hence using the formula 



i Hik -qk){Pk-dk) 
{Qk - If 

-duqk + {-lYqu^'Pk 



k=l i=l 



(Qk - iy+' 



2m-l 



^.^2m-l ^ ^ /'2m l^|^.Q,^2m-l-, 
w — n V J 



j=0 



and grouping in powers of /3 we have 



+ 



Ti 



^2m 2m-2 



+ 



E 



(2m)! (2m -1)! 
(2m-l)!/32™-i-i 



/3 



AT+l 

2m-l Vifc 



(gf^^ -qk){Pk-dk) 



{qk - ly 



(2m- 1)! ^ j!(2m-l-j)! 

2jn, m,— 1 2m— 1 



7 + 1 2^ 2^ ' 

•' k=l i=l 



{Qk - 1; 



j+i 



+ 



+ 



K' 



(2m- 



1)1 E E 



-duqu + (-l)^gr+>fc 



ir\2m—l 



ilk - i; 



j+i 



A*0 



(43) 



k=l i=l 

Now consider T2. Using binomial formula and equalities (5) the expression for T2 we reduce 
to the form 



T2 = E^W- 
7=0 



2(2m - 1) 



13 



1^' {h/3f"'-^-^{-iy 



z 



E 



{h/3) 



m—j f 1 \ m+j — 1 



1)' 



AT 



-r,Y.^b]br^'-'- (44) 



2^0' + 1)K2^ - 1 - j)! 2 ^ (m + J - l)!(m - j)! ^ 
And for the right hand side fm[P] of equation (4) (see (10)) using binomial formula we get 



UP] 



{h(5) 



m— 1 



(2m) 



2m j^ "i-i 



(/i/3)2™-i--'(-l)J+i 
(j + l)!(2m-l-j)! 



E 



(/i/3)™-J'(-l)"'^-' 



(45) 



Then substituting (43), (44) into (40) taking into account (45) after some simplifications for 
the difference /„ [/?] — g [P] we get 



fmm-9[(i] 



l^2m+ 



(2m - 1) 



o2m-l V^ 



™-l t^N+1 



ilk -Qk)idk-Pk) 



2m-2 



l^2m -™ -(2m-l)!/32™-i-^' 
(2m -1)! ^ j!(2m-l-j)! 



fe=i 

m— 1 J 



(gfc - 1) 



^^■+1 + 'v" y^ -4gfc + (-i)'gf+Vfc ^^Q, 



j + 



fc=l i=l 



(gfc - 1)^ 



+1 



+ 



1 ^ {hp)"'-^{-l)'^+^ 

2 ^ (m + j - l)!(m - j)! 



Af 



m + j 



Y.^bWir-''-' 



7=0 



^2m ™-l 2m-l ^_^ , ^ ^ ^^^Af+i„ 



(2m - 1) 



_1^1 Z^ Z^ 



-4gfc + {-lyqr'Pk 



ir\2m—l 



A*0 



k=l i=l 

Hence note that /m[/3] — g[P] is polynomial of degree (2m — 1) of [/3] = /i/3, i.e. 

2m-l 



i=o 



Here 



/l2 



m-1 , ]v+ 



(2m~l) 



bi + 



1)\ ^ 



fc=l 



fe-i)2 



(-ly 



J ' 2j\(2m-l-j)\ \ 2m-j 



N 



ECbM 



2m-j-l 



7=0 



m—1 2m— 1 



fe^'" "y^' 'y^' -t^fcgfc+(-i)'gf"^>fc Atn2m-i 

"(2m-l)! ^ ^ 



fc=l J=l 



fe-l)*+^ 



: /\in2m—l i 



1 



2(2m-l)! I 2m 



N 



ECbM 



2m-l 



7=0 



for j = 2m — 1, 

for m < j < 2m — 2, 

for 1 < j < m — 1, 

for j = 0, 



(46) 
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where 



l^2m-j 



B. 



2m— j 



j\(2m — j — \)\ \ 2m — j 



2m— J — 1 m— 1 

;+ E E 






^ig2m-jr-l 



i=l k=l 

On the other hand from (6) for difference fm[P] — g[(3] the following is true 



fmm - gm = Pm-im. 



If 



m-l j 

EE 

fc=l i=l 



4gfc + i-iy^'qr'Pk 



{qk - 1; 



j+i 



A^O^' 



-^, J = l,m-1, 
J + 1 



y^^ (9f^^-?fc)(4-Pfc) 



0, 



(47) 

(48) 
(49) 



k=l 



{qk - ir 

then equality (47) is take placed. From (46) and (47) we find unknown polynomial Pm-i[f3] in 
the system (6)- (8) 

m— 1 

^ 1^' (50) 



Pm-im = J2aA/3y 



Thus, for unknowns dk and pk {k = l,m — 1) we have got system of linear equations (48) -(49). 



Later, from (8) for a = 2, m — 1, using equalities (37), (39), (25), (26) and (49) after some 
simplifications we obtain 



'"-I " / , i . / , y+l AT+l^ J ^Ar+ 



h'*' E E 



k=l i=0 



d,4 + (~ir\r'pk 



+E 

i=2 



aW 



B 



A*0"- 



4gr^ + (-l)^+^gfcPfc 
(1 - qkY^' 



A*0" + 



(j -!)!(« + 1-j)! 



J 



J sp sp "fcgfc +1--'-; qkPk ^jQj-i 



^-qk 



\i+l 



0. 



(51) 



k=l i=0 

The left hand side of the equation (51) is the polynomial of degree (a + 1) with respect to h. 
Then from (51) obviously that all coefficients of this polynomial are zero, i.e. 



k=l 1=0 



+1 



[^-qk 



\i+l 



j + r 



j = l,a-l, 



(52) 



J^ J2 dkql + {-iy^'qr'Pk ,.^. ^J^f: ^'^^^ + {-^y^'1^P^.^ 



k=l i=0 



[i-qk 



k=l i=0 



;i - qkY^' 



A*0°, a = 2, m-l. (53) 
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Thus for unknowns dk and pk we have got the system (52), (53). 

Using lemmas 2, 3 for the left hand side of the system (53) we obtain 

fc=l j=0 ^ ' k=l 1=0 ^ ' 

where a = 2, 3, ..., m — 1. 

Hence taking into account (48) we get 

g^M±H^:^^&A-0<' = (-ir«|^, „ = 2,3.....™-l. (54) 

fc=l i=0 ^ ^^> " ^ 

Since for a = 2£, £ = 1,2, ... Bernoulli numbers -Bq,+i = -B2£+i are zero, then combining (54) 
and (53) we have 

yydj^^^H-}r^^.^.^B^ a = 2,3,...,m-l. (55) 

ttt^ (1-g.)^-^^ « + l' 

Comparing systems (52) and (55) it is easy to see, that from (52) for j = 1 we get only one 
equation with respect to unknowns rf^ and p^. And for j = 2, ..., a — 1 obtained equations from 
system (52) are part of the system (55). 

Thus, from (48), (52), (55) we get the following system 



k=\ i=\ 



EE ^'^t;- J^ "^ A-O^-^, .-1,2,...,^-!, (56) 



-'^^grvir^+g.i^...,,.^ 



y y -m, y-^j j ,kPk ^., ^ ^^ j = 1,2, ...,m - 1. (57) 



k=l j=o 

Subtracting equality (57) from (56) gives 

y:V (^^-^^^-^Ci^'^-^^^ A-0^- = 0, , = l,2,...,m-l. (58) 

From system (58) for j = 1 we get equation (49). Therefore for finding unknowns d^ and pk it 
is sufficiently to solve the system (56), (57). 

Further we investigate the system (58). After some transformations homogenous system 
(58) we reduce to the form 

E '""^%'::n"'--"' -°. .-u.....™-i. (59) 



(* - 1)' 
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where dk — pk, k = 1,2, ..., m — 1 are unknowns. Obviously, that the main matrix of the system 
(59) is the Vandermonde type matrix, i.e. system (59) has a unique solution, which identically 
zero. This means 

dk=Pk, k = l,2,...,m-l. (60) 

Then keeping in mind (60) from (37), (38), (39) and (56) we respectively get the following 

system 



k=l 



C[0] = h 



C[N] = h 



i- -y d,^-^^^^\ 



Qk + (-ij Qk Mr,j _ ^j + l 



ydk y,''\.:' A-0^' = ^^, j = l,2,...,m-l. 
Theorem 3 is proved. 

References 

[1] Sobolev S.L. Introduction to the theory of cubature formulas. Moscow, Nauka , 1974, 808p. 

[2] Shadimetov Kh.M. Discrete analogue of the operator ^^ and its construction // 
Quary of computational and applied mathematics. - Tashkent, 1985. - pp. 22-35. 
ArXiv.l001.0556vl[NA.Math]. 

[3] Hamming R.V. Numerical methods. Moscow, Nauka, 1968, 400 p. 

[4] Gelfond A.O. Calculus of finiute difference. -Moscow, Nauka, 1967. -367 p. 

[5] Frobenius. Uber Bernoullische Zahlen und Eulersche Polynomen Sitzungsberichte der 
Preussischen Akademic de Wissenschaften. 1910. s.809. 



17 



